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ABSTRACT 

Recent Spitzer observations of the globular cluster M15 detected dust associated with its intracluster 
medium. Surprisingly, these observations imply that the dust must be very short-lived compared to 
the time since the last Galactic plane crossing of the cluster. Here we propose a simple mechanism to 
explain this short lifetime. We argue that the kinetic energy of the material ejected during a stellar 
collision may be sufficient to remove the gas and dust entirely from a cluster, or to remove the gas as 
a wind, in addition to partially destroying the dust. By calculating the rate of stellar collisions using 
an iV-body model for the cluster, we find remarkable agreement between the average time between 
collisions and the inferred dust lifetime in this cluster, suggesting a possible close relation between 
the two phenomena. Furthermore, we also obtain the birthrate of blue stragglers formed through 
collisions in M15. By comparing with the observed number of blue stragglers, we derive an upper 
limit for their average lifetime which turns out to be consistent with recent model calculations, thereby 
lending further support to our model. 

Subject headings: stellar dynamics - shock waves - blue stragglers - stars: mass loss- globular clusters: 
general - globular clusters: individual (M15) 



1. INTRODUCTION 

In recent years, infrared space telescopes like IRAS, 
ISO and especially Spitzer made it possible to investigate 
not only the stellar content of globular clusters (GCs) but 
also their intracluster medium (ICM), consisting of gas 
and dust lost by giants. Evans et al. (2003) were the first 
to clearly detect thermal dust emission from the core of 
the GC M15 (NGC 7078) with ISO. They found a dust 
mass an order of magnitude lower than one would expect 
based on the number, mass-loss rates and lifetimes of 
horizontal branch (HB) stars, even at the extremely low 
metallicity of M15. Subsequent, more sensitive measure- 
ments with Spitzer resulted in a dust mass of 9 x 10~ 4 M 
indicating an accumulation time-span not much longer 
than 10 6 yr (Boyer et al. 2006) . While such a short time 
is commonly explained by ram-pressure stripping during 
Galactic plane passages, the time since the last passage 
for M15 is about 4 x 10 7 yr, an order of magnitude larger. 
Therefore, additional processes have been suggested to 
explain short ICM lifetimes, including blow-out of nova 
explosions, fast winds from stars, relativistic winds from 
millisecond pulsars (Spergel 1991), and radiative ejection 
by the strong radiation field in M15 (Smith 1999). 

Another possibility that has not yet been investigated 
is related to the outcome of stellar collisions. From hy- 
drodynamic simulations it is known that colliding main- 
sequence stars (MSS) lose typically ~ 1% of their mass 
(see, e.g., Lombardi et al. 1996). As this material is 
ejected with speeds typically a few times the escape 
speed from the stellar surface, it will transfer a significant 
amount of energy to the ICM and possibly reduce the 
amount of observable dust considerably. For instance, if 
the gas, released through a collision between two 0.5 M 
MSSs, has a mass of 0.02 M and leaves with twice the 
escape velocity (« 1200 km s" 1 ), it possesses enough ki- 
netic energy to accelerate more than 15 M Q , the esti- 
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mated total mass of ICM gas in M15 (Boyer et al. 2006), 
up to a speed of 40kms~ , the escape speed from the 
center of M15. 

The inferred dust lifetime in the ICM would then 
require that such a stellar collision happens every ~ 
1 Myr. Given that M15 has a central density of n > 
10 6 M Q pc~ 3 and a velocity dispersion of a « lOkrns -1 
(see, e.g., McNamara et al. 2004), the collision time 

T coU = 7xl0 11 yr (^f 1 ) ( TooKT^) for a 1 M o MSS 

is ~ 10 10 yr (Binney & Tremaine 1987). This results in 
at least one collision every < 5 Myr for the more than 
2000 stars in the core (e.g., Dull et al. 1997). From this 
rough estimate one can already see that this mechanism 
is indeed promising. 

However, there are several factors in M15 that compli- 
cate a better estimate of the collision time. Firstly, in 
GCs as old as M15, all MSSs have masses below 0.8 M Q , 
while initially more massive stars have produced rem- 
nants like white dwarfs (WDs) and neutron stars (NSs), 
which will have different interaction rates. Secondly, 
because of mass segregation the lower-mass MSSs may 
concentrate in lower density regions, while more massive 
WDs and NSs concentrate near the center. Finally, it 
is not a priori clear that enough energy from the gas 
ejected in a stellar collision can be effectively transferred 
to the ICM since losses through radiative shocks may be 
important. 

We address the issue of mass-segregation considering a 
specific model of M15 given in McNamara et al. (2004), 
which is based on an iV-body model of Baumgardt et 
al. (2003). This allows us to estimate collision times in- 
volving different stellar populations in the cluster. Fur- 
thermore, we obtain a minimum formation rate for blue 
stragglers (BSs) through MS-MS mergers. Using the 
number of observed BSs in the core of M15 we estimate 
an average BS lifetime and compare it with evolution- 
ary models, as an extra check on our basic model (f2"]l. 
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In Sj3] we discuss how the ejected material from a stel- 
lar collision interacts with the ICM and how it can clear 
signatures of dust emission. We conclude in 

2. COLLISION TIMES 

We define a collision between two stars with radii i?i 
and i? 2 to occur whenever their distance d < R\ + R 2 . 

(1 2) 

The average local collision time TV,', (r) for one star of 
type "1" to collide with a star of type "2" can be written 
as (compare to Binney & Tremaine 1987 their eq. 8-122) 
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where ni.2 is the number density of field stars "2" and 
stars of type "1", Mi^ their masses, and a the one- 
dimensional velocity dispersion (assuming a Maxwellian 
velocity distribution for both species with er(r) = 
\J (ci(^) 2 + &2{r) 2 )/2) . To get the total number of col- 
lisions per unit time in the cluster between these two 
species, eq. Q]is integrated over the whole cluster, 



dNcoiijot 
dt 



= 47T 



I 

drr 2 m(r) - — — 
T co u{r) 



(2) 



where r is the radial position in the cluster. 

In order to account for a continuous stellar mass spec- 
trum and the mass-radius relationship of stars, we take 
local averages Ri,%(r) and Mi j2 (r) at position r. We find 
that in the N-body model for M15 these profiles, as well 
as 71^2 (r), can be well represented by power-laws over 
a sufficiently large range of 0.025 pc < r < 2pc. Inside 
of 0.025 pc there are almost no MSSs, while outside of 
« 1 pc the contribution of the integrand in eq. [2] becomes 
rapidly negligible. We do not consider collisions between 
MSSs and NSs as the NS retention fraction in GCs is 
expected to be low (Pfahl & Rappaport 2001; Dull et al. 
1997). We also do not consider collisions between giants 
and MSSs as the escape velocity at the surface of a gi- 
ant, and even more so the expected energy of the ejected 
material, is an order of magnitude lower than for MSSs. 
For simplicity we choose a constant a = llkms -1 , as 
a does not vary much within 1 pc (±lkms _ ; compare 
Dull et al. 1997) and also agrees with the value obtained 
by McNamara & Baumgardt (2004) for M15 for a similar 
region (r < 0.3' » 0.8pc). The slight variations within 
1 pc are accounted for in the error estimates of the col- 
lision rates. With all quantities given as power-laws, we 
solve eq. [2] analytically. Our calculations are based on 
the results of iV-body simulations by Baumgardt et al. 
(2003), which were scaled to fit the velocity dispersion 
profile of M15 in McNamara et al. (2004). Their model 
consisted of initially 130, 072 stars with a realistic mass 
spectrum, and included a treatment of stellar evolution 
and the Galactic tidal field. They also take into account 
velocity kicks imparted to NSs at birth and consider two 
extreme cases, one where all NSs are retained, and one 
where all NSs are removed from the cluster. Here we 
only consider the latter case, since, as mentioned before, 
the actual NS retention fraction is expected to be very 
low. 

In Table [1] the fit parameters for n, R, R 2 and M 



TABLE 1 

Parameters of the power- law fits for the mass (M), 

DENSITY (n), AND STELLAR RADIUS (R) FOR WDS AND MSSS. THE 
POWER-LAWS ARE OF THE FORM a = do {r/ro) b WITH Tq = 1 pc. 





a 


b 


n-MS 


(3.9 ±0.2) x 10 3 pc- 3 


-1.620 ±0.006 


n WD 


(3.3 ±0.1) x 10 3 pc- 3 


-2.280 ±0.009 


Rms 


(0.51 ± 0.05) R Q 


-0.143 ± 0.008 


MS 


(0.32 ± 0.05) Rq 


-0.26 ± 0.01 


M MS 


(0.49 ± 0.02) M 


-0.105 ± 0.002 


M WD 


(0.74 ± 0.02) M 


-0.117 ±0.002 



for WDs and MSSs in M15 are shown. As can be seen, 
through mass segregation, MSSs are more abundant out- 
side 1 pc while the more massive WDs dominate the cen- 
tral 0.75 pc where collisions are most likely to happen. 
As a consequence, we should expect more collisions be- 
tween WDs and MSSs than between MSSs. From eq. O 
we obtain dN co u^MS-WD/dt = (2 ± 1) x 10 _7 yr^ 1 and 
dN coU , MS -Ms/dt = (4±l)x lO^yr -1 , where N coUtX -Y 
is the number of collisions between X and Y. It follows 
that, given the current state of M15, we expect one colli- 
sion every (5 ± 2) x 10 6 yr between a WD and a MSS and 
every (3 ± 1) x 10 7 yr a collision between two MSSs, re- 
sulting in one collision every (4 ±2) x 10 6 yr that releases 
high-velocity gas into the ICM. 

It is rather remarkable how closely this timescale coin- 
cides with the estimated ICM dust lifetime, suggesting a 
direct connection between dust clearing and stellar colli- 
sions. 

Using a similar approach, we estimate the BS forma- 
tion rate, defining BSs as two merged MSSs with a com- 
bined mass > 1 M Q , significantly larger than the MS 
turn-off mass (0.8 M ). After binning all MSSs into 
0.1 Mq bins, we determine the collision rate between all 
mass bins with combined mass > 1 Mq. This also allows 
us to determine the expected mass spectrum for BSs for 
comparison with future observations. We obtain a total 
collision rate of (6.4 ±0.7) x 10 _3 Myr _1 which, together 
with the 6 — 7 BSs observed in M15, implies an average 
BS lifetime of about 1 Gyr. This is also consistent with 
recent BS evolution models (e.g., Leigh et al. 2007; Sills 
et al. 2001). However, we note that this is certainly an 
upper limit, given that the presence of binaries would 
increase the formation rate through binary mergers and 
resonant interactions (Fregeau et al. 2004; Mapelli et al. 
2004). For example, for a Plummer sphere with a Kroupa 
mass function ranging from 0.3 to 3.0 Mq and a binary 
fraction of 30%, the rate of collisions mediated by binary- 
single and binary-binary interactions can be ~ 1 order 
of magnitude larger than from single-single interactions 
(Chatterjee et al. 2008). In Figure Q] the mass distribu- 
tion for BSs is shown. As expected it decreases for higher 
mass BSs, since the number density for the lower mass 
MSSs is much higher compared to that for the close to 
turn-off MSSs. 

3. INTERACTION OF COLLISION EJECTA WITH THE ICM 

Interactions between very fast moving gas and a low- 
density interstellar ambient medium have been exten- 
sively investigated in the context of supernova rem- 
nants. Their evolution progresses in several distinct 
stages (Chevalier 1977): the "ejecta dominated" (ED) 
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Fig. 1. — Estimated mass distribution of BS candidates created 
via collisions between two MSSs with combined mass > 0.8Mq. 
rii, b, and N are the number of BSs in the ith bin, the bin size, 
and the total number of BSs, respectively. 



stage, characterized by a freely expanding blast-shock 
wave until the mass of the swept up material is com- 
parable to the mass of the ejecta, the "Scdov- Taylor" 
(ST) stage, where the blast-wave expands adiabatically, 
and the "pressure-driven snowplow" (PDS) stage, where 
a thin shell forms which "snowplows" through the am- 
bient medium, driven by the pressure of the hot in- 
terior in addition to its own momentum (Cox 1972; 
Cioffi et al. 1988). Cioffi et al. (1988) obtained a sim- 
ple offset power-law solution that describes the kinemat- 
ics of the blast wave in the PDS stage, which can be 

written as v s = vpp>s {Rs/ Rpds) where v s and 
R s are the velocity and radial position of the shock 
front, respectively, relative to the site of the collision, 
and vpjjs and Rpds are the velocity and radial po- 
sition of the shock front at the transition from the 
ST stage to the PDS stage. The transition values are 



given by v PDS = 413n ' Q 



1/7^3/14^1/14 



51 



kms 1 and Rpds 



14.0 E 5 [ ( m n 1 pc, where £51 is the kinetic energy 
of the ejected material in units of 10 51 ergs, no the den- 
sity of the ambient medium in cm -3 , and £ m is the heavy 
element abundance relative to solar abundances. An in- 
teresting property of the flow of the post-shock gas is 
that in the PDS stage its mean velocity has the same 
value as the shock speed (Cioffi et al. 1988). For our 
problem this means that, provided v s is larger than the 
cluster escape speed after the shock swept up all of the 
ICM, the shocked gas can entirely escape the cluster. 

It is now interesting to see if such an ICM removal 
mechanism might be applicable to M15. For this esti- 
mate we assume an ICM mass of 15 M Q , which is ex- 
pected given the amount of dust detected by Boyer et 
al. (2006) and the low metallicity of M15 (Cm = 10" 24 ). 
In order to estimate E 51 for collisions between MSS we 
can use the results of Lombardi et al. (2002) (their Table 
3), obtaining £51 = 1 — 11 x 10" 4 . Unfortunately, there 
are no similar results for MS-WD collisions. However, 
given that the mass loss is presumably between 15% to 



50% (Ruffert & Mueller 1991; Ruffert 1992) and that 
the velocity of the ejecta should be of the order of the 
escape speed of the MSS (« 500 kms" 1 ) we obtain a 
similar range for E 5 i. Note that, although the ejecta 
are ejected nearly isotropically (Lombardi et al. 1996), 
the ICM is likely to have a more irregular structure, as 
the patchy dust emission in M15 indicates. Thus, only a 
fraction of the released energy might actually be trans- 
ferred to the ICM. As the 3D structure of the ICM is 
unknown, we limit our analysis to a fiducial ejecta en- 
ergy of -E51 = 3 x 10" 4 but also determine the minimum 
E51 value to remove all of the ICM gas from the cluster, 
noting that these values are to be understood as "ef- 
fective" energies ramming into the ICM gas. Similarly, 
estimates for the ICM density are also rather uncertain. 
If we assume that all the gas is contained within a radius 
of 1 — 2 pc (the approximate radial position of the dust 
emission), the density of a 15 M Q gas cloud would be be- 
tween 20— 150 cm -3 . For simplicity, we also assume that 
the shock, once it leaves the ST stage, remains radiative, 
which might not necessarily be the case for the low tem- 
peratures (T < 10 5 K) and shock speeds (< 50kms _1 ) 
when the shock has swept up most of the ICM, as cool- 
ing is less efficient in this regime (Sutherland & Dopita 
1993). 

Using E 51 = 3 x 10" 4 and n = 20 - 150, cm" 3 we 
obtain v s w 13 — 15 kms" 1 after the shock has swept 
up 15 M Q . Since this velocity is lower than the escape 
speed of ~ 40 kms" 1 for M15, it follows that the energy 
of the ejected material from one MSS collision may not 
be sufficient to remove the ICM completely. In fact, only 
for E51 > 8 x 10" 4 , which is at the upper end of the 
estimated ejecta energy interval, the shocked ICM gas 
attains a velocity sufficiently large to leave the cluster. 

On the other hand, the shock also heats up the gas to 
high temperatures, (« 10 — 15 x 10 4 K at those shock 
speeds). Gas at such temperatures expands beyond a 
critical radius and flows out of the cluster as a wind, 
reducing the amounts of ionized gas down to less than 
1 M© (Knapp et al. 1996). Assuming that the dust fol- 
lows the gas, the dust would therefore leave the cluster on 
a timescale as short as ~ 10 5 yr. In fact, the wind should 
be even stronger in our scenario than for the static model 
considered in Knapp et al. (1996), since here the gas itself 
has a considerable outwards speed through the shock. 

So far, we assumed the existence of rather large 
amounts of gas, based on the amount of observed dust 
and M15's extremely low metallicity. However, searches 
for gas in M15 have had very limited success, and Smith 
et al. (1995) estimate an upper limit for the total ICM 
mass of about 3 M Q . The reason for the much larger ob- 
served dust-to-gas ratio is not well understood (see, e.g., 
van Loon et al. 2006). As the dust-to-gas ratio of the 
material lost in the winds of red giants should scale with 
the metallicity of the stars (van Loon et al. 2005), thus 
resulting in more than 10 M of gas, it appears that addi- 
tional processes may be at work that remove the gas more 
easily than the dust. Nevertheless, it is also possible that 
most of the ICM is in molecular form as the CO-to-H2 
conversion factor is not known for such low metallici- 
ties and extreme radiative environments (van Loon et al. 
2006). Repeating the previous calculation for an ICM 
mass of M = 3 M Q and n = 3 — 30 cm" 3 accordingly, we 
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obtain v s — 52 — 63kms~ 1 for E$i = 3x 10~ 4 , while for 
£51 = 2 x 1CT 4 v s = 39 — 47kms~ 1 which is close to and 
larger than the cluster escape speed of 40km s _1 . In this 
case, the ejecta of one MSS collision would likely be able 
to accelerate this gas out of the cluster. 

As a final caveat, we note that it is not very clear 
whether the dust will follow the rather low-density ion- 
ized gas. This strongly depends, among other quantities, 
on the dust grain properties and their electric poten- 
tial relative to the ionized gas in a rather complicated 
way (Draine & Salpeter 1979). For instance, Nozawa 
et al. (2006) and similarly Slavin et al. (2004) found in 
their simulations of shocks driven into dusty interstellar 
medium that, while small grains with sizes of w 0.01 /im 
get destroyed by sputtering and grain-grain collisions, 
only grains with sizes w 0.1 /jm are actually dragged 
along with the gas, while grains with sizes > 1 /iin remain 
almost unaffected and do not follow the shock wave. On 
the other hand, if the size distribution of the dust grains 
is similar to the one for the local interstellar medium, 
we see (e.g., Mathis 1996) that most of the dust mass is 
in grains with sizes of w 0.1 /xm. This means that, even 
if not all dust particles follow the gas flow, we can nev- 
ertheless expect that most of the dust mass will remain 
sufficiently well coupled to the gas and consequently be 
removed. 

4. CONCLUSIONS 

In this Letter we proposed a new mechanism to explain 
the relatively short lifetimes of the ICM dust in a dense 
GC, developing our arguments in detail for the case of 
M15. By calculating the rate of stellar collisions using the 
detailed model for M15 by McNamara et al. (2004), we 
find a remarkable, close agreement between the average 
time between collisions and the inferred dust lifetime of 
~ 10 6 yr (Boyer et al. 2006) in this cluster, pointing to 
a direct link between the two phenomena. We argue 
that the kinetic energy of the material ejected during 
a stellar collision may be sufficient to remove the dust 
from the cluster, depending on the assumed ICM mass, 
either directly by accelerating dust and gas to velocities 



larger than the cluster escape speed, or indirectly, by 
accelerating and heating the gas, which then expands 
and leaves the cluster as a wind, carrying the dust along 
with it. Although there are some uncertainties as to how 
well the dust will couple to the gas, especially at the 
low shock speeds expected for this problem, there are 
some indications from simulations of supernova remnants 
that might support sufficient coupling (e.g., Nozawa et 
al. 2006). In addition, at least some dust grains can 
also be efficiently destroyed by grain-grain collisions or 
sputtering (Slavin et al. 2004; Nozawa et al. 2006), which 
further helps to reduce the amount of observable dust 
in the cluster. On the other hand, the results of these 
studies may not be directly applicable to our scenario 
because, e.g., the intense UV field present in a cluster 
like M15 could strongly affect the electric potential of 
the grains and, therefore, their coupling to the ionized 
gas (Draine & Salpeter 1979). 

With a detailed model for M15 we were also able to cal- 
culate the formation rate and mass distribution of BSs 
through MS-MS collisions. By comparing with the ob- 
served number of BSs in the cluster, we derive an upper 
limit to their average lifetime of ~ 1 Gyr, consistent with 
current stellar structure and evolution models for BSs 
(e.g., Leigh et al. 2007). 

We conclude that the interaction of ejected gas from 
stellar collisions with the ICM will strongly influence the 
observable signal of the dust in the ICM, and, given the 
remarkable agreement between dust lifetimes and colli- 
sion times in M15, this represents a promising mechanism 
to explain the very short dust and ICM lifetimes in GCs. 
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